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Abstract 

We analyze the phase diagram of a system of spin-1/2 Heisenberg antiferromagnetic chains 
interacting through a zig-zag coupling, also called zig-zag ladders. Using bosonization 
techniques we study how a spin-gap or more generally plateaux in magnetization curves 
arise in different situations. While for coupled XX Z chains, one has to deal with a recently 
discovered chiral perturbation, the coupling term which is present for normal ladders is 
restored by an external magnetic field, dimerization or the presence of charge carriers. 
We then proceed with a numerical investigation of the phase diagram of two coupled 
Heisenberg chains in the presence of a magnetic field. Unusual behaviour is found for 
ferromagnetic coupled antiferromagnetic chains. Finally, for three (and more) legs one 
can choose different inequivalent types of coupling between the chains. We find that the 
three-leg ladder can exhibit a spin-gap and/or non-trivial plateaux in the magnetization 
curve whose appearance strongly depends on the choice of coupling. 
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I. INTRODUCTION 

In the last few years, the study of quasi-one dimensional magnets has become intense. 
One of the main reasons is the appearance of real materials which can be well approximated 
by one-dimensional models. An important class corresponds to the so-called spin ladder 
materials, such as the compounds Sri_ x Cu x 02 and La4 + 4a;Cu8+2a:0 u+sa which are closely 
related to high-T c compounds (for reviews see e.g. |l|||) or the organic two- leg ladder material 
Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 (see e.g. §). 

In addition, the similarities in the normal-state properties of high-temperature super- 
conductors and ladder cuprates make these latter systems valuable laboratories from both 
the theoretical and experimental point of view. 

Due to the low dimensionality, quantum fluctuations are crucial and because of this, these 
systems exhibit a variety of interesting phenomena such as the appearance of plateaux in 
magnetization curves, an issue that has received a lot of attention recently (see e.g. [P|-pD|]). 

The study of spin ladder systems with different topologies of couplings has been also 
motivated both from the experimental and the theoretical side. In particular, the zig-zag 
coupling between quantum spin chains has received much theoretical attention [pl|-p5|. In 



the case of antiferromagnetic (AF) couplings, the zig-zag array introduces frustration which 
makes the study of these systems much more difficult. Apart from being a possible approach 
to study the two-dimensional triangular lattice, this topology of couplings is realized in a 
number of quasi-one dimensional compounds, such as Cs 2 CuCl 4 PB[, KCuCl 3 and TlCuCl 3 



27j as well as NH4CUCI3 [28|. A (two-dimensional) zig-zag arrangement is also present in 
SrCu0 2 . At room temperature, the materials KCUCI3, TICUCI3 and NH4CUCI3 are isostruc- 
tural and can be described by an alternating two-leg zig-zag ladder [E7|,E3]. Experimentally, 



one observes a zero magnetization plateau in the low-temperature magnetization process of 
KCUCI3 and TICUCI3 0], whereas for NH4CUCI3 plateaux in the magnetization curve are 
observed at 1/4 and 3/4 of the saturation magnetization p8| . While the former is in good 



agreement with the theoretical predictions (which will be discussed in detail in this paper), 
the explanation of the latter is still unclear (even though explanations have already been 
proposed [[£J). A detailed analysis of related structures such as those studied in the present 
paper can be expected to be also useful for interpreting these experiments. 

Zig-zag coupled chains also have several special properties which render them an inter- 
esting problem from a purely theoretical point of view. Firstly, unlike many other systems 
they do not have a simple 'strong-coupling' limit where the system is decoupled into finite 
clusters of spins. If there is such a decoupling limit, the possible values of the magnetization 
are clearly quantized in this limit and using perturbation theory one can easily understand 
the appearance of magnetization plateaux in the strong-coupling region [|TI],[T^,[nj . For ex- 
ample, for the usual spin ladders, one finds a quantization condition on the magnetization 
(M) (which we normalize to saturation value 1) for the appearance of a plateau in the 
magnetization curve 0,0 (compare also |6|,|^,|l^, |il||ilfl ): 



SV{l-{M))eZ. (1.1) 

Here S is the spin on each site and V is the volume of a translationally invariant unit cell. 
For iV-leg ladders, V = IN where / is the period of explicit or spontaneous breaking of 
translational symmetry in the magnetized groundstate. 



Since zig-zag coupled chains with no or weak dimerization do not have a simple strong- 
coupling limit, it is not immediately clear if the condition for the appearance of plateaux in 
zig-zag ladders will also be given by ( |1 . 1| ) . The fact that the zig-zag coupling is frustrating 



is a further reason why the quantization condition on the magnetization or e.g. the univer- 
sality classes of the transitions at the plateau boundaries might be different. From a field 
theoretical point of view, zig-zag coupling is also interesting because at zero field it cancels 
the most relevant coupling term for the usual ladders and instead one has to deal with a 
chirally asymmetric perturbation |24|| . 

After mentioning all these possibilities for a different behaviour we should, however, 
immediately point out that one of our main conclusions is going to be that zig-zag coupling 
is not very different from the ordinary one. In particular, the most relevant interaction 
term is recovered in many situations like the presence of a magnetic field, dimerization or 
charge carriers. We will also find that all observed magnetization plateaux have a natural 
interpretation in terms of the quantization condition ( |1 . 1| ) . 

Motivation for studying these systems arises also from other fields. For example, spin 



ladders arise in the study of gated Josephson junction arrays [pOfl . Even more, Ref. |22 
pointed out a close analogy between the two-leg zig-zag ladder and the doped Kondo lat- 
tice model. This model in turn is believed to be relevant to the phenomenon of colossal 
magnetoresistance in metallic oxides and has therefore received renewed attention recently 
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The main focus of the present paper are N coupled Heisenberg chains with a dimerized 
zig-zag coupling (see also Fig. [I]) 

N L N L 

H (N) = J Y / Y / s i , x - s i>x+1 + </'££ 4* ■ ((i + 5 ) s i+ i,x + (i - s) swm+i) - hJ2 s* >x . 

i=l x=l i=l x=l i,x 

(1.2) 

For compactness of presentation we have written scalar products here, but we will also 
consider the case where an XX Z anisotropy A is introduced in the obvious way. We will 
furthermore discuss a similar system including dimerization along the chains and also charge 
degrees of freedom, i.e. Hubbard zig-zag ladders. In eq. ( |1.2| ) there is also some ambiguity 
in writing the interchain coupling, in particular in combination with specifying boundary 
conditions for N > 2 - an issue to which we shall return later. 

The two-leg zig-zag ladder can be reinterpreted as a Heisenberg chain with next-nearest- 
neighbour interactions. In this interpretation, the study of this system has a long history 



going back at least to |j32[ . One intriguing observation regarding this system is that for 
special parameters the groundstate takes a very simple form and correlation functions can 
be computed exactly [0,@ • 



This paper is organized as follows: In Section II we examine the Hamiltonian ( |1.2|) at 
5 = h = without and with XX Z anisotropy, using non-Abelian and Abelian bosonization, 
respectively. Section III is devoted to a field-theoretical analysis of the effect of various 
modifications: An external magnetic field, dimerization or doping with charge carriers. In 
Section IV we then numerically investigate the magnetization process of the two-leg lad- 
der ( |1.2| ) and discuss the various commensurate and incommensurate phases appearing for 
J' > as well as J' < 0. In Section V we investigate the transition to saturation of the 



two-leg zig-zag ladder in some detail. Finally, in Section VI we numerically compute mag- 
netization curves for several variants of the three-leg zig-zag ladder. An appendix contains 
some supplementary material. 
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FIG. 1. Generic structure of the ladders considered in the present paper. 



II. FIELD THEORY FOR ZERO MAGNETIC FIELD 



A. The SU{2) symmetric case 

First, we reexamine the simplest case of the SU(2) symmetric two-leg zig-zag ladder that 
has been studied previously in |22||23| . 

We start by reviewing some of the earlier results. In the weak (interchain) coupling 
limit each of the chains can be described by the level 1 Wess-Zumino-Witten (WZW) model 



35|,|3"6|j with the action given by 
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where the superscript i labels the different chains, g l takes values in the Lie group 577(2) 
and V[g] is the Wess-Zumino term 
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(2.2) 



with Y a three-dimensional manifold having M 2 as boundary. 
The bosonized expression for the spin operator is given by 

S\x) = T R + Jl + B{-lfti{ag l ) 



(2.3) 



where J% a = -^tr{{d + g l ){g i )- 1 cx a ), JJ a = ^tr((^ i )- 1 9_^V a ) are the WZW currents satis- 
fying a su(2) 1 Kac-Moody algebra and B is a non-universal constant. 



The marginally irrelevant perturbation term 

-X^fdxJ^-Jl, A>0, (2.4) 

i 

which is responsible for the logarithmic corrections in the correlators in the case of a single 
spin chain is usually discarded when there are more relevant terms (such as those arising in 
the normal ladders from the interchain coupling). In the present case, though, it should be 
kept since the zig-zag interchain coupling gives rise only to marginal perturbations. 
More precisely, the zig-zag coupling 

Hint = J 2^ S x ■ (S x + S x+1 ), (2.5) 

X 

leads to the (marginal) current- current interaction 

H mt = aj dx(Jl + Ji) • (4 + f[), (2.6) 

where a is small and proportional to the interchain coupling J'. It is positive for anti- 
ferromagnetic coupling and negative for a ferromagnetic coupling. Discarding non-Lorentz 
invariant terms, which do not contribute to the one-loop renormalization group (RG) equa- 
tions, one can rewrite the Hamiltonian in the following way 

H = H 1 +H 2 -PJ2J d < T R ■■%)+«/ dxJE ■ J[ (2.7) 



t=i 



where (3 = X + a and J 1 = J 1 + J 2 which satisfies a level 2 Kac- Moody algebra. Equation 
( p.7| ) is the Hamiltonian considered in p2]j23| . By a one-loop RG analysis one finds an 
exponentially small gap for antiferromagnetic zig-zag coupling and a massless regime for a 
ferromagnetic coupling. However, in [^31 it was shown that there is another (marginal) term 

7 J dx (triag 1 ) ■ d x tr{°9 2 ) ~ tr(<?# 2 ) • ^tr^ 1 )) , (2.8) 

with 7 oc J', that could change the behaviour of the system in the ferromagnetic regime. 



Now we proceed to discuss how the RG equations of |23| are changed when the chirally 
asymmetric perturbation term (|2.8j ) is included in the non-Abelian bosonization framework. 
Such a computation was actually already performed in |3J] using a fermionic representation. 
Still, our result will be complementary to |24j] and is specially appropriate to analyze the 



SU(2) symmetric case and generalize it to iV > 2 since this symmetry is explicit in the 
present analysis. 

It is easy to see that another SU(2) invariant operator given by 

7 ' / dx (tr^dMg 2 ) - tr^dMo 1 )) , (2-9) 

has to be included as a counterterm in the action. Using the operator product expansion of 
these operators and a rescaling of the coupling constants, one obtains the RG equations: 
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a 2 + 7 2 , 

-f3 2 + 2a(3 + 7 2 , 

2 7 a + 2 7 /3 + 2 7 'a , 

67a - &y'/3 + 67a . (2.10) 

The RG equations obtained in [E3| can be recovered from here by just putting 7 = 7' = 0. 
It can be shown that for an AF coupling the presence of these 7 and 7' terms does not affect 
the qualitative behaviour of the flow. The system is driven to a strong coupling regime, 



corresponding to the massive behaviour mentioned in ||23| . The way terms like (|2.8|) affect 



the nature of the massive excitations in this strong coupling regime is nevertheless not so 
clear. 

These equations can be generalized to an arbitrary number of coupled chains, provided 
all the chains are coupled in the same way. This applies in particular to iV = 3 with peri- 
odic boundary conditions (PBC). One can then conclude that three chains with a periodic 
arrangement of antiferromagnetic coupling J' > are massive at least in the weak-coupling 
region J' <C J. Due to the dimensions of the perturbing operators involved, this kind of 
coupling can be considered as weaker than the normal coupling, in the sense that the growth 
of the associated coupling is governed by a smaller exponent in the present case. This shows 
the importance of a careful treatment of the marginal interactions. This should help to 
clarify the problem of the presence of a gap in the weak-coupling region for the usual N = 3 
ladder with PBC. 

The situation for a ferromagnetic coupling is more subtle. One can see that the trivial 
massless infrared fixed point a = (3 = is stable only if 7 = 7' = 0. Then the presence 
of the twist term affects the large scale behaviour of the system, preventing it to reach the 
trivial massless point. So, the one-loop RG flows for weak J' < and J' > become similar. 
This has already been noticed in [E3| where it was argued that the dimerized phase seems 
to extend into the ferromagnetic region. One should however keep in mind that this is just 
a one-loop calculation and non-perturbative effects could change the large scale behaviour 
of the system. 

We conclude the present section stressing that all these results are easily generalizable 
to the weak-coupling regime of an arbitrary number N of coupled chains, provided that in 
the continuum limit all chains are coupled in an equal manner. It should be noted, however, 
that on the lattice, inequivalent versions of completely symmetric zig-zag interchain coupling 
exist for iV > 2. After taking the continuum limit, these differences manifest themselves in 
different signs for interaction terms of the type ( |2.8| ), i.e. most of the originally completely 
symmetric boundary conditions are not symmetric anymore after taking the continuum limit. 
Note also that the cases of PBC where not all pair-couplings are present (which is the case 
for N > 3) or open boundary conditions (OBC) for N > 2 are much more subtle and it is 
not clear if the above results apply also to them. 



B. The XXZ case 

Let us consider now the addition of an XXZ anisotropy. In this case, the SU(2) symme- 
try of the system is broken and it is then more appropriate to adopt the Abelian bosonization 



approach. Below we follow the notations of |36| , |22[| in order to simplify making the connec 



tion with the non- Abelian description in terms of WZW fields of the previous section (this 
implies some minor changes of conventions such as a rescaling of the fields with respect to 
our study of the ordinary spin ladders |I3| ). 

In the Abelian bosonization language, each chain is described by a compactified free 
bosonic field 0* with its dynamics governed by [] 

H= \S dx ( vK ^^ 2 + ^(<^) 2 ) ( 2 - n ) 

The field <p % and its dual 4> l are given by the sum and difference of the lightcone components, 
respectively. The constant K governs the conformal dimensions of the bosonic vertex oper- 
ators and can be obtained exactly from the Bethe ansatz solution of the XXZ chain (see 



e.g. [13] for a detailed summary). We have K = 1 for the SU(2) symmetric case (A = 1) 
and it is related to the radius R of [13j by K^ 1 = 2ttR 2 . 
In terms of these fields, the spin operators read 

Sf x = -^d x (t) 1 + a : cos(2^x + v 7 ^) : +^r 1 , (2.12) 

v27r 2 



Sg. = (-1)- : e^^ 1 (b C0S (2k F x + V2^(f) 1 ) + c) : , (2.13) 

where the colons denote normal ordering with respect to the groundstate with magnetization 
(M l ). The Fermi momentum k F is related to the magnetization of the ith chain as k l F = 
(1 — (M 1 ))tt/2. The effect of an XXZ anisotropy and/or the external magnetic field (to be 
discussed later) is then to modify the scaling dimensions of the physical fields through K and 
the commensurability properties of the spin operators, as can be seen from ( [2. 12 ), ( p. 13 ). 



The non-universal constants a, b and c can be in general computed numerically (see e.g. 
|37j, for the case of zero magnetic field) and in particular the constant b has been obtained 



exactly in 

At zero magnetization, (i.e. k F = tt/2), the first and second terms in each of these 
equations correspond to the components of the non-oscillatory (J l R + J l L ) and oscillatory 
(tr(ag 1 )) terms in ( |2.3|) . 

Let us consider now a two-leg zig-zag XXZ ladder. The perturbation terms in the case 
of zero magnetization can be separated into three classes: 

i) Terms quadratic in the derivatives 

ad x <p\x)d x <p\x) } (2.14) 



1 Note that a factor (4i?) 2 is missing in front of the II 2 term in eqs. (2.2) and (3.4) of 13] 



that can be absorbed into a renormalization of the compactification radii, once we 
diagonalize the derivative part going to the new variables ^ = (0 1 ± 2 )/v / 2- The K 
parameters are then renormalized as 

K± = K(l T 2J'K/(Jn) + 0((J'/J) 2 )). (2.15) 

ii) The other contributions from current-current intrachain and interchain interactions 
can be rewritten as 



— A cos(v47T0 + ) cos(v47T0 ) 

+ a cos(v / 47T0 _ ) (cos(v / 47r0 + ) + cos(\/47r</r)) . ( 2 -16) 

The dimensions of these operators are given by 

K + + K^ and -L + K+, (2.17) 

respectively (the third operator in ( |2.16| ) is always irrelevant). The A term (which is 
the standard current- current term for the individual chains) is irrelevant for A < 1, 
while the a term is irrelevant for a ferromagnetic interchain coupling, except for K — 1 
(and J' = 0) which corresponds to the SU(2) symmetric point where it is marginal. 
For an AF coupling, the a term is relevant for A > 1. 



iii) We also have chirally asymmetric terms |23| (those corresponding to ( |2.8| )) which now 
read 



-d x (p + (x)sm(V^(p~); +d x (p- (x) sin(v^7T0 + ) ; -d x ^j) + {x) sin(v / 4^0") . (2.18) 
Their dimensions are 

1 + K_, 1 + K+, and 1 + -^- . (2.19) 

If we consider a small XX Z anisotropy with A -C 1, both K± increase and the term 

d x 4> + (x) sm(y/47T(f)^) 
will now become the most relevant for sufficiently small J'. 



Following arguments similar to the ones of [|24j it follows that at most the symmetric 
{4> + ) sector is massless. For A > 1, also this field acquires a mass, but stays massless 
for A C 1. The above analysis cannot directly be applied to the region A « 1, because 
higher loops should be included. Moreover, at the point A = 1, the SU(2) symmetry is not 
directly explicit in this treatment. However, according to the one- loop RG analysis in f2H}| 
and Section II. A of the present paper, a gap seems to open at A = 1 for small \J'\ ^ 0, 
regardless of its sign. 



III. STABILITY OF THE ZIG-ZAG COUPLING AND MAGNETIZATION 

PLATEAUX 

In the present section we will study different mechanisms due to which the zig-zag inter- 
chain coupling becomes unstable against the perpendicular interchain coupling. 

The outcome is that only under very special circumstances the zig-zag coupling is stable, 
while under the action of an external magnetic field, the addition of rung dimerization or 
the presence of charge carriers, an effective perpendicular coupling is generated. 

A. XXZ in a magnetic field 

When we add an external magnetic field to the XXZ zig-zag ladder, the situation is 
more subtle, since the staggered terms do not necessarily cancel. Indeed, we show in this 
section that new non-oscillating commensurate terms arise. 

The current-current terms become 



sin/ci?cos(v47T0 ) (sin kp cos(v4tt(J) + + Akpx) + sin(/c j p — V47T0 )) , (3.1) 

and the other terms become 

- d x (j) + (x) sin(2fcp - v 7 ^^) -= (d x 4> 2 (x) m\(Ak F x + 2k F + v / 47T0 + ) 

V 2 v 

-d x <f>\x) sin(4fc F x - 2k F + v / 4^0 + )) - d x 4> + (x) sin(v / 4^^) . (3.2) 



The main difference with the case at (M) = is that now we have a very relevant term 
which is proportional to 



cos k F cos(4k F x + V4it(J) + ) + cos A^cos^i? — V47r</> ) (3.3) 



with dimensions K + and K_. Note that the first term in (3J3) will disappear since it is 
incommensurate for (M) ^ 0. Hence this term gives a mass to the <ft~ field, leaving only the 
symmetric field massless. 

We can then make the following statement about the phase diagram of two zig-zag 
coupled XXZ spin chains with A < 1: For non-zero magnetization, and in all cases, one of 
the degrees of freedom is massive and the other massless, leaving a c = 1 theory. 



With some modifications to the argumentation of ||13| , this can be easily generalized 
to N weakly coupled zig-zag chains, provided all the chains are coupled together: The 
generalization of the most relevant interaction term (|3.3|) to the case of iV-leg ladders is: 



y^ j Ai cos 2 k F : cos(Axk F + V27r(0i + 4>j)) ■ +X2 cos kp : cos(kp — v27r(0j — 0,,)) : \ . (3.4) 



1 -j 



As for the case N = 2, the coupling constants Aj essentially correspond to the coupling J' 
between the chains: Aj ~ J'/ J, but have a non-trivial dependence on (M): Aj — > for 
(M) — > 0. The Gaussian part of the Hamiltonian is now given by: 



H {N) = J^\lZ{ K {dJ l (x)f + I (X0*(*)) 2 } + ^ £ {dMx)) (d x ^(x)) 



N 

1 / ' 1 ,■ -. ■ ( y\ 

(3.5) 

—, 1 v. i\ > i\ ■ ■ 

L i=l v ' 1,] 



where A ~ J' /J. As for the case of the normal coupled ladder, the last term produces a 
shift of the compactification radii of the fields and plays a crucial role in the opening of 
non-trivial plateaux. In arriving to the Hamiltonian ( |3.5| ) we have discarded a constant 
term and absorbed a term linear in the derivatives of the free bosons into a redefinition of 
the applied magnetic field. One has also to include the generalization to N chains of the 
current-current term and the twist term. While obtaining the expression for the first one is 
an easy task, writing down the twist term for generic N has some subtleties. Keeping only 
commensurate terms, a naive generalization of ( |2.18|) would be: 

]T -d x (fa + fa) sin(2A; F - V^ifa - fa)) - d x {fa + fa) sm(v^r(& - fa)) • (3-6) 

i>j 



If we take into account only the most relevant perturbation term (|3.4j) , the analysis is similar 
to the one of the normal ladder given in [0|. In particular, one can radiatively generate 
"iV-Umklapp" terms like cos [2xJ2iLi kp + V^nNipD), where i/jd = h= YliLi fa, which can 
give rise to the appearance of a plateau if ( |1.1[ ) is satisfied with V = N, S = 1/2. 



The zero-loop analysis for the opening of such plateaux is basically the same as in |]T3 
Just the coefficient of the interaction arising from the smooth part acquires an extra factor 
of two. Via the Gaussian part of the Hamiltonian this replaces J' by 2 J' in the dimension 
formula for the iV-Umklapp term. In this way, we obtain for example, at A = 1, J' c « 0.045 J 
for the (M) = 1/3 plateau at N = 3 and J' c « 0.35J for (M) = 1/2 at iV = 4 and also 
for (M) = 1/5 at iV = 5. At this level one could expect the magnetization curves of 
zig-zag ladders to be very similar to the ones of normal ladders. There are however some 
subtleties which are not captured by this analysis. There are inequivalent ways of coupling a 
number N > 3 chains in a periodic zig-zag way. These inequivalent ways of PBC couplings 
correspond in the weak coupling field theoretical model to changing the relative signs, or 
what is equivalent, permuting chain indices in the expression of the twist term (|3.(j|). Since 



the zero-loop treatment presented above does not take into account the twist term, it is 
obvious that we have to consider loop corrections to take into account this effect. A detailed 
RG treatment of this model taking into account the current- current and twist term is a 
difficult task beyond the scope of this paper. However, exact diagonalization of finite chains 
for a strong enough J' /J ratio will show that these inequivalent couplings can give rise to 
different behaviour in the presence of a magnetic field. Another problem is the extension of 
these results to OBC where we encounter the same limitations as for the normal ladders. We 



refer the reader to [|T3[ for a detailed discussion of this point. We just mention here that the 
results above cannot be directly extrapolated to configurations with OBC. In this case one 
should therefore complement the present field theoretical analysis by other methods such as 
exact diagonalization. 



B. Dimerization in XXZ ladders 

Another way to generate an effective perpendicular coupling is to include dimerization 
in the zig-zag coupling. Let us consider the Hamiltonian in (|1.2| ) with non-zero rung dimer- 
ization 5. In the bosonized language, this gives rise to the perturbation term 

H dim = J'EEft- 1 )^ • ((! +5)(-l) x fi i+ i + (l-8)(-iy +1 n i+1 



10 



2J ' S J2J2™i '™i+l> 



(3.7) 



where rij is the staggered component of the spin operator at zero magnetization [n\ 
cos(\/27r0*), nf = exp(±i\/27r0 1 ))- 

This perturbation term can be rewritten as 



H dim = A J2 J2 [~ cos (V^r(0j + 0i+i)) + cos (v / 47r(0 



?i+i 



+2 cos (v / 47r(0i - 



^*+i. 



(3.8) 



where A oc J' 5. 

The effective model is then very similar to the one for a ladder with a normal perpen- 
dicular coupling. In particular, using the formulae of [|F|], one can study the opening of 
plateaux in the magnetization curve as a function of S and J'. The N = 2 version of the 
dimerized zig-zag ladder has been analyzed in detail in |12j and the magnetization plateau 
with (M) = 1/2 predicted there was also observed numerically [39|J^0[ . 

Another way to add dimerization is along the legs of the chains as in Fig. |2|, which gives 
rise to additional terms. They read in the presence of a magnetic field 



v 



5 J2 J2(~ 1 ) X cos{2k F x + V2^(j) 



(3.9) 



i=l x 



which is incommensurate at non-zero magnetization, but can contribute to radiatively gen- 
erated commensurate terms as we discuss below (see also (|20f). 



J (1-5) J(l+<5) 




FIG. 2. Dimerization along the legs. 
Taking as an example N = 2, one can show that the operator 



(3.10) 



is radiatively generated for non-zero magnetization, and it is commensurate for (M) = 1/2. 
This operator arises from a combined effect of the interchain coupling and the dimerization 
along the chains and hence A' in (|3.10| ) is proportional to J', 5 and (M). This operator 
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has dimension K and is hence relevant, implying the existence of a plateau at 1/2 of the 
saturation magnetization. 

A related analysis was performed in [BJ but in the regime where the interchain coupling 
J 1 — > — oo, which cannot be studied within the perturbative approach used here. In fact, 
in that case the existence of the plateau was associated with a different operator of higher 
dimension, which could become relevant only under certain conditions. 

In the present case, the existence of the plateau is independent of the values of the 
microscopic parameters, since the dimension of the operator ( |3.10| ) is always smaller than 
2. Furthermore, for a two-leg zig-zag ladder there is only one way to introduce dimerization 



along the legs and therefore in contrast to the ordinary ladders [f20[ cancellations do not 
occur. 

An amusing observation which may also be of some relevance for experiments is that 
one recovers the two-leg zig-zag ladder with dimerized interchain coupling as the first-order 
effective Hamiltonian for the two-leg ladder with dimerization along the chains (see Fig. 0) 
in the limit of 8 — > 1. Details are given in an appendix. The main point is that dimerized 
interchain coupling and dimerization along the legs break translational symmetry in different 
ways: For the former case a translationally invariant unit cell of the Hamiltonian contains 
two spins while in the latter case it contains four. In view of the quantization condition 



[Tp, this suggests that dimerization along the legs can give rise to more plateaux than if 
only the interchain coupling is dimerized. 

A similar analysis as the one of this section can be performed for generic N. 

C. Doping with charge carriers 

Let us consider now a generalization of the system studied in Section II by adding charge 
degrees of freedom. We start by describing the Hamiltonian of interacting electrons in one 



dimension [36,4] 



H = ~17j2( C l+l,» C x,a + H.c.) + Uj2 c t,+ c x,+ct,-c x ,- (3.11) 

" x,a x 

where c^ and c are electron creation and annihilation operators and a — ±. For positive 
U and at half filling, the charge sector is massive and the spin sector for large U can be 



described by the Heisenberg Hamiltonian ||36|| . Here, we will analyze the large scale behaviour 
of two identical systems coupled in zig-zag, away from half filling. Since we are now keeping 
the SU(2) symmetry of the spin sector (neither magnetic field nor XX Z anisotropy will be 
considered), the continuum limit of the Hamiltonian at non-zero doping can be written as: 

H= \! dx (^c(<9.0 c ) 2 + -^(a^ c ) 2 ) + H wzw {g s ) -\jdxJ R ■ J L , (3.12) 

where the first term of the Hamiltonian stands for the charge sector and the WZW term 
describes the spin sector. The charge and spin density operators are given respectively by: 

p{x) = jn + Jl + const. sm(2kpx + V2n(j) c ) tr(g) + const. cos(Akpx + v8vr0 c ) (3.13) 

and 
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— * — * — * 1 

S{x) = Jr + Jl + const. sm(2k F x + v 27r0 c ) tr(ag) , (3-14) 

where j and J are the U(l) and SU(2) currents of the charge and spin sector respectively. 
The Fermi momentum kp is now related to the charge sector and we keep our spin sector at 
zero magnetization. It is known that K c — 1 for U — and K c = 1/2 for U = oo ||42|| . 

We again consider two copies of the system and study perturbatively the zig-zag coupling 
between them. More precisely, we consider a coupling term which involves both charge and 
spin densities 



tfi* = E (jp'M + pI+x) + J'slisl + S^) 



The presence of a further direct hopping term between the chains would require a more 
detailed analysis. Indeed, first the Gaussian part must be diagonalized including this direct 
hopping term and then, for example, one can treat perturbatively terms like H int . This 
is beyond the scope of the present article where we just concentrate on showing that the 



presence of a weak coupling between Hubbard chains like in |41j generalized to the zig-zag 
configuration can give rise to results similar to those observed for normal ladders. 

In the continuum description, this interaction term can be written in two different pieces. 
We first have the current-current interaction term given by 

y(ii + Jl)U 2 R + Jl) + J'(Jh + Jl) ■ (4 + Jl) (3-15) 

The U(l) current- current term can be completely absorbed in the Gaussian charge Hamil- 
tonian by a rescaling of v c and K c 

where the indices ± stand for the symmetric and antisymmetric fields (<p l c + (jfy/y2~ and 
((f>l — (fi?)/y/2 respectively. The SU(2) current term is identical to the one of Section II. A 
and we know then that it plays a role only for positive J'. 

There is however another commensurate perturbation present away from half filling which 
is given by 



coskp cos(v47T0 c ) [U' const. tr((7 1 ) tr(g 2 ) + J' const, tr^g 1 ) ■ tr(ag 2 )) . (3-17) 



In the present case, the chirally asymmetric terms are irrelevant in the RG sense. 

In the weak interchain coupling limit we are considering, U' is much smaller than U. 
Then U' /(2ttv~) < 1 — K 2 and this term is relevant giving a mass to all the spin degrees 
of freedom as well as to the antisymmetric charge field <fi~. Thus, only the field 0+ remains 
massless. 

We see from (|3.17|) that also in this case an effective perpendicular coupling is generated, 
signaling the instability of the zig-zag interchain coupling also against charge doping. 

Since the perturbation terms in ( 3 . 1 7|) are the same as those that appear in the normal 



coupling studied in [H], it is then natural to expect that the IR behaviour of the charge 



density wave and superconducting correlation functions will be the same in the present case. 
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IV. NUMERICAL ANALYSIS OF THE TWO-LEG CASE 

In this section we consider the N = 2 version of ( |1.2| ) without dimerization, i.e. 5 = 0. 
It is then useful to think of the two-leg zig-zag ladder as a single chain with next-nearest 
neighbour interaction. So, for N = 2 the Hamiltonian (|1.2| ) can be recast in the form 
(writing explicitly an XX Z anisotropy A) 

H = J Z^ \ ^x^x+1 + 9 \^x S x+ i + S x S X+1 J > 

x=l L Z J 

+ ^ Z^ j ^^x^x+2 + 9 V x x + 2 ~*~ x x+2j r (4-1) 



L 

x=l 



Here L denotes the total volume of the system. In this section we will always assume J > 0. 
To avoid frustration also in the limit J' — > where we find two weakly coupled chains, we 
choose L to be a multiple of 4. 

In the formulation (|4.1|) essentially all spatial symmetries are manifestly implemented 
by a one-site translation x — > x + 1. So, one can use Fourier transforms to simplify the 
determination of the spectrum. Since the magnetic field h is coupled to a conserved order 
parameter in fl4.1|) , we can relate all quantities at a field h to those at h — - the results to 
be reported below are all obtained from computations with h = 0. 

There is already a number of exact diagonalization studies for the two-leg XX Z zig-zag 
ladder (or equivalently the Heisenberg chain with next- nearest-neighbour interactions) in 
a magnetic field P5|-^5|,p5|] . However, there are still some regions in the parameter space 
and aspects which have not been studied in great detail, such that further numerical inves- 
tigations seem worthwhile. For a numerical analysis we concentrate on the isotropic point 
A = 1. It turns out that groundstate momenta can be incommensurate in the presence of 
the magnetic field. Such a feature is interesting in its own right and reminds us of recent ob- 
servations made for the 5 = 1 chain with biquadratic interaction |50| , |51"f , but also gives rise 
to technical complications. This incommensurability is one reason why the computations to 
be reported below needed substantially more CPU time than analogous computations for 



conventional ladders |I(J , even though we used an improved version of the program employed 
loc.cit. 

In order to scan the whole range of coupling constants we choose a normalization such 
that J + | J'\ = 1. J' = then corresponds to two decoupled chains, J' = 1 to a single 
antiferromagnetic chain and J' = —1 to a single ferromagnetic chain. The resulting mag- 
netic phase diagram is shown in Fig. |3|. Here, the lines show the magnetic fields where 
the magnetization jumps between two different values that are realized at a given system 
size. The conclusions are schematically summarized in the inset of Fig. |3|. Regions with 
a ferromagnetically polarized groundstate are denoted by an 'F'. The other regions will be 
explained in the following discussion of our results. 

Let us look first at the antiferromagnetic region J' > 0. Here spins flip one after the 
other. The groundstate momenta are in general incommensurate (i.e. not multiples of 7r/2) 
in the region IC1 (0 < J' < 4 J), though for the small lattice sizes accessible to us, this 
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FIG. 3. Magnetic phase diagram with the choice of normalization J + \J'\ = 1. The lines are 
for L = 24 (full), L = 20 (dashed). L = 16 (dotted), L = 12 (long dashed-dotted) and L = 8 
(dashed-dotted). The inset shows a schematic version which is discussed in the text. 



does not show up in the region of small J'. For h = 0, a study of the static structure factor 
exhibited a transition to incommensurate behaviour at J' /J ~ 1.92075 [|52[| . 

For h > 0, the onset of incommensurability in the groundstate momenta was determined 
in j|7|], and |x| determined a transition in the groundstate using a different criterion. In 
the incommensurate region IC1 of Fig. |3|, the lines are irregular which is mainly due to the 
lattice discretization of the momenta. In the region CI, i.e. for J' > 4 J, all groundstate 
momenta are commensurate and convergence with system size is good. 

A gap [21],p^1 can be anticipated in Fig. |3| for J < J' < 2 J, but apart from that there is no 
evidence for any non-trivial plateaux. In fact, non-trivial plateaux have not been observed 
at the various points which have been studied over the past decade |^,[£|[|4],[|(],(|j|] . These 
observations can be nicely interpreted in terms of the quantization condition ( |1 . 1| ) . In view 
of the mapping to a single chain ([4.1|), one should substitute iV = 1 in (|1.1|). The gap then 
arises by spontaneous breaking of this enhanced translational symmetry to / = 2. Non- 
trivial plateaux would then require / > 2, which at least in the two- leg zig-zag ladder does 
not seem to be permitted. 

The magnetization process on the ferromagnetic side has already been looked at some 



time ago fr45| , though no definite conclusion was reached due to problems of resolution. 
Indeed, the ferromagnetic side is quite a bit different from the antiferromagnetic one, even 
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though plateaux do neither seem to occur here. In the region C2 (—3 J/2 < J' < 0), spins flip 
in pairs and the groundstate momenta of those states that participate in the magnetization 
process are commensurate (actually or ±7r). The states with an odd number of spins 
pointing up (or down) can have incommensurate momenta, but they do not participate in 
the magnetization process. 

At J' = —4 J a transition to a completely polarized ferromagnet takes place. The precise 
location of the transition point can be attributed to the fact that at J' = —4 J an exact S = 



groundstate, the uniformly distributed RVB state [M j can be written down and is degenerate 
with the ferromagnetic state. The intermediate region IC2 (i.e. roughly — 4 J < J' < —3 J/2) 
is rather complicated: Here the number of spins flipping simultaneously at a given system 
size changes with J'/ J. Furthermore, groundstates with incommensurate momenta (i.e. 
momenta not an integer multiple of n) do participate in the finite-size magnetization process. 
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FIG. 4. Magnetization curves at J' = — 3 J with lines for L = 24 (full), L = 20 (dashed). 
L = 16 (dotted), L = 12 (long dashed-dotted) and L = 8 (dashed-dotted). The bold line is a 
sketch of the form expected in the thermodynamic limit. 



The magnetization process in the region IC2 is illustrated in Fig. |] by the finite-size 
magnetization curves at J' = —3J. Here, the number of spins that flip simultaneously in 
a finite-size system varies between 1 and 3 (the number of spins flipping simultaneously 
increases as one approaches J' = —4 J). Strong non-monotonic finite-size effects can be 
observed in particular at the smaller system sizes. Nevertheless, a reasonable approximation 
to the limit L — >• oo seems to be obtained applying the procedure of ]55[ to L = 20 and 
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L = 24, i.e. by connecting the midpoints of the steps in the finite-size magnetization curves. 
This yields the bold line in Fig. |]. The behaviour at small fields is somewhat speculative 
since the finite-size gaps are non-monotonic such that it is not possible to extrapolate them. 
We have therefore simply assumed a vanishing gap in the limit L — > oo. 

The region close to saturation in Fig. |] is on comparably safe grounds. In fact, we have 
data for up to L = 48 which has been taken into account for the bold curve. In the case 
J' = —3 J, this data yields no evidence for more than three spins flipping simultaneously at 
the transition (M) — > 1. Therefore we are confident that the magnetization curve becomes 
smooth in the thermodynamic limit (though very steep at the transition (M) — > 1). In 
particular, we think that our intermediate phase IC2 is different from the metamagnetic 



one observed in [||| in a different parameter region, even though at first sight they bear 



some resemblance. The fact that we are not aware of any evidence for incommensurate 



momenta arising in the metamagnetic phase ||8[ also suggests that there are two distinct 



phases preceding a transition to a ferromagnet in different parameter regions. 

We would like to conclude the present discussion with a few remarks on the spin- gap. 
Although it has already been investigated with high accuracy for J' > applying the density 



matrix renormalization group to large system sizes pl~| , |22"| , |2"5|| , there are still some interest- 
ing problems. Firstly, only two of these investigations f2"T| , |2"2"|] deal with the incommensurate 
region and there the results differ by more than 10%. Secondly, no such investigations 
were performed for J' < 0. In these two regions, 5* = 1 excitations (corresponding to the 
lowest lines in Fig. ^| for J' > 0) have an incommensurate momentum at the minimum of 
their dispersion. Since the system size limits the resolution in momentum-space, one has 
additional non-monotonic finite-size effects. We have tried to overcome this with an inter- 
polation using a Fourier transform of the dispersion. However, even/odd momentum effects 
(which are characteristic for scattering states) obscure the true incommensurate minimum 
and preclude such an analysis. According to preliminary investigations, such an approach 
is applicable to S — 1/2 excitations (the fundamental spinon-type excitations) not only in 
the commensurate region |25| but also in the incommensurate region. 

The question of a spin-gap is particularly interesting for J' < 0. In Section II. A we have 
confirmed the conclusion of |23| that a one-loop RG analysis suggests a non-zero spin-gap for 
small J' < 0. However, for J' < there is no region where a gap in the magnetic excitations 
is obvious from our numerical data. The data for the S z = 2 excitations (corresponding to 
the lowest lines for J' < in Fig. |3|) actually fits well to the form 

E(L)~~. (4.2) 

Such a form is characteristic for a conformally invariant and thus gapless situation (see 
e.g. [|5B|). For this reason we have not drawn a gap in the schematic inset of Fig. |3| for 
J' < 0. Nevertheless, it remains a challenge to find numerical evidence for a spin-gap in 
some magnetic excitations and look for possible massless sectors which should be organized 
into su(2) 1 representations. 

V. THE TRANSITION TO SATURATION IN THE TWO-LEG LADDER 
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A. The antiferromagnetic side 

Now let us look at the transition to saturation in the lattice-version of the N = 2-leg 
zig-zag ladder. 

First we recall the computation of the transition field h nc [Q. One can immediately 
write down the energy of a fully polarized (ferromagnetic) state (with h = 0) 

E aat . = jA ( J' + J) . (5.1) 

Using a Fourier transform, also the excitation energy of a single spin-wave above this fully 
polarized state is readily computed as 

g u (p) = -A{J' + J) + J'cosp+Jcos{2p). (5.2) 

To determine the critical magnetic field h nc associated to the transition (M) — ► 1, we have 
to minimize (|5.2|) . One finds the minimum at 

p min =h-«»- 1 (£) forO<J'<4J, (53) 

' 1 vr for J' > 4 J. 

An immediate consequence of this simple result is that the groundstate of the S z = L — 
1 sector has an incommensurate momentum for < J' < 4J. This is presumably the 



simplest example of how the two competing interactions in ( f4.1|) lead to the phase with 
incommensurate groundstate momenta which we discussed in the previous section. 
Insertion of ( |5.3|) into (|5.2j ) directly leads to the upper critical field 



h JAJ'+(A+1)JH-^ for0<J'<4J, (54) 

\(A+1)J'+(A-1)J forJ'>4J. V ; 

Now we will go beyond EHf and discuss the nature of the transition (M) — *■ 1. To this 



end, it is useful to look at two-spinwave excitations |57j • After a Fourier transformation, one 
arrives at a matrix problem in the distance of the two flipped spins. We omit the explicit 
expressions for the matrices that we have used to study the behaviour of the magnetization 
curve close to saturation. This is equivalent to studying the finite-size behaviour of the 
two-spinwave groundstate energy £ 2s at h — h uc , since 1 — (M) = 4/L and £ 2s — h — h uc . 

For the two-leg zig-zag ladder, the eigenvalue problem in the two-spinwave subspace 
can be interpreted as a fourth-order difference equation with suitable boundary conditions 
and thus can in principle be solved by a further Fourier transformation. For our purposes 
it turned out to be sufficient to fix the center-of-mass momentum p and then perform a 
numerical diagonalization for system sizes up to L ~ 150. We have looked at A = 1 
and the following values of the coupling constants: J' = 2 J, J' = 4 J, J' = 4.1J and 
J' = 4cos(7r/12) J. These values were selected since for them ( |5.3|) can be realized exactly 
for suitable choices of the system size L and one can thus avoid further non-monotonic 
finite-size effects which could arise from the lattice discretization in momentum space. In all 
cases that we have studied, the minimum energy of the two-spinwave excitation was found 
in the sector with center-of-mass momentum p = 0. 



The numerical diagonalization determines the critical field associated to the transition 
from two flipped spins to one flipped spin. For J' ^ 4 J, this critical field is compatible with 
the universal DN-PT behaviour |58| , |59| 

/i uc - /i ~ (1 - (M» 2 . (5.5) 

Just at J' = 4J we find a different behaviour which is much better described by [] 

h uc -h~(l-{M))*. (5.6) 



To understand why just the point J' = 4J leads to a different exponent, it is instructive 
to look at the behaviour of £\ s near the minimum p m i n . . From (|5.2|) and ( |5.3|) we find 

[ (4J - j y j+J V + (9(x 3 ) forJ'<4J, 
£i s Gw + x) - £ u {pmin.) = \i^ + 0{x Q ) for J 1 = 4J, (5.7) 

[ (f - 2 j) x 2 + £>(x 4 ) for J' > 4J. 

We see that the dispersion is quadratic everywhere except for the point J' = 4J where it is 
quartic. These exponents in the one-spinwave dispersion are in one-to-one correspondence 
with those in (|5.5|) and ( |5.6j ), respectively. So, the behaviour of the magnetization for 
(M) — > 1 can be explained by the band structure. Note that ( |5.7| ) is independent of A, and 
so should be the exponents corresponding to ( |5.5|) and (|5.6j ), respectively (as long as A > 0). 

B. The ferromagnetic side 

Now we turn to the transition to saturation on the ferromagnetic side J' < 0. Since we 
know from the numerical investigation that spins flip in pairs (as long as \J'\ is not too large), 
we have to solve a non-trivial matrix problem already to determine the critical field h uc in 
the ferromagnetic regime J' < 0. For \J'\ < 2.5 J, the minimum of the dispersion of the 
two-spinon state is now located at p = it. In the region J' < — 2.66J, the minimum moves 
away form p = n, i.e. it becomes incommensurate. In the following we restrict ourselves 
to the region of sufficiently small \J'\ such that spins flip in pairs and only commensurate 
groundstates participate in the magnetization process for (M) — > 1. 

Some values of transition fields h uc are listed in Table |. These have actually been 
obtained on chains with a few hundred sites. Nevertheless, all given digits should be those 
of the thermodynamic limit. 

To check the asymptotic behaviour of the magnetization curve, one further needs at least 
the transition field ft.4^2 at which two spins flip to reduce the number of spins deviating 
from the ferromagnetic state from four to two. Table | lists transition fields /i4^2 computed 
numerically on chains with 36 < L < 72. Since 1 — (M) ~ 1/L, we can then test for the 
universal behaviour (|5.5| ) by forming the expression 



2 The exponent 4 has also been observed in [46| at J' / J = 4. 
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L\ h n , - h 



4^2 



(5.8) 



which should converge to a constant if ( |5.5j ) is satisfied. Indeed, the expression (|5.8|) appears 
to converge to a constant with incresing L for all cases listed in Table |. So, we find no 
counterevidence against the universal DN-PT behaviour on the ferromagnetic side either, 
although we have no compelling argument in its favour. 



J'/J 





-3/13 


-1/3 


-5/11 


-3/5 


-7/9 


-1 


-9/7 


-5/3 


-11/5 


L 










/l4_> 


2/J 










36 


1.96595 


1.77368 


1.69888 


1.60799 


1.50460 


1.38523 


1.24387 


1.07095 


0.85064 


0.55455 


40 


1.97272 


1.78215 


1.70116 


1.60981 


1.50630 


1.38682 


1.24527 


1.07208 


0.85145 


0.55495 


44 


1.97766 


1.78427 


1.70263 


1.61109 


1.50753 


1.38796 


1.24624 


1.07287 


0.85201 


0.55522 


48 


1.98137 


1.78568 


1.70366 


1.61204 


1.50844 


1.38879 


1.24695 


1.07343 


0.85241 


0.55542 


52 


1.98423 


1.78666 


1.70441 


1.61277 


1.50914 


1.38941 


1.24748 


1.07385 


0.85271 


0.55556 


56 


1.98648 


1.78736 


1.70499 


1.61335 


1.50968 


1.38988 


1.24788 


1.07418 


0.85294 


0.55567 


60 


1.98828 


1.78788 


1.70546 


1.61381 


1.51010 


1.39025 


1.24819 


1.07443 


0.85312 


0.55576 


64 


1.98974 


1.78828 


1.70583 


1.61418 


1.51044 


1.39054 


1.24844 


1.07463 


0.85326 


0.55582 


68 


1.99094 


1.78859 


1.70615 


1.61448 


1.51071 


1.39078 


1.24864 


1.07479 


0.85338 


0.55588 


72 


1.99195 


1.78885 


1.70641 


1.61473 


1.51093 


1.39097 


1.24880 


1.07492 


0.85347 


0.55592 


h uc /J 


2 


1.79086 


1.70833 


1.61648 


1.51250 


1.39236 


1.25000 


1.07589 


0.85417 


0.55625 



TABLE I. Values of h A ^ 2 at 
ferromagnetic state takes place at 
of h uc which were estimated using 



which a transition from four to two 
a given size L. The last row contains 
much larger systems. 



spins deviating 
auxiliary data, i 



from the 
. e. values 



We expect that the transition (M) — > 1 remains second order if we increase \J'\ further 
into the regime where more than two spins flip at the same time (this is in contrast to the 
metamagnetic transition studied in |Q). At least in this region the very steep behaviour 
of the magnetization curve (see Fig. ^) could mean a breakdown of the DN-PT universal 
behaviour ( |5.5| ). However, due to the flipping of several spins at the same time, we are not 
able to investigate the transition to saturation in detail numerically in this region. 



VI. THE THREE-LEG CASE 

Now we proceed with a discussion of A^ = 3-leg zig-zag ladders on the lattice, though 
in less detail than for the two-leg ladder. For A^ = 3 there is a large number of possible 
boundary conditions for the coupling between the chains. In particular for periodic boundary 
conditions (PBC), there are several possibilities already for AT = 3, of which none is naturally 
singled out in the case of zig-zag coupling. We will discuss three types of PBC in addition to 
open boundary conditions along the rungs. To elucidate these different types of boundary 
conditions it may be useful to consider the case J' = J where the three-leg zig-zag ladder can 
be considered as a strip of the triangular lattice Heisenberg antiferromagnet. The triangular 
lattice has three sublattices and it is possible to consider boundary conditions that do or 
do not respect this sublattice structure. We will call the PBC that respect this sublattice 
structure at J' = J 'PBC of type A', while those that identify different sublattices 'PBC of 
type B and C. 
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The iV = 3-leg zig-zag ladder with PBC of type B and C can be rewritten as a single 
Heisenberg chain with interactions up to distances of three, respective four sites: 



( 1 1 

H = J Z^ \ ^n^n+i + 9 \S n S n+1 + S n S n+1 \ + AS n S n+r + - yS n S n+r + S n S n+r 

3L ^i n 3L 

+ ^ Z^ 1 ^^n^n+3 + 9 \pn ^n+3 + ^n ^n+3j [ ~ " Z^ ^n ' (^-l) 

n.=l L Z J n=l 

Here the second interaction term goes over a distance r = 2 for PBC of type B, while PBC 
of type C are characterized by r = 4. 



Due to (|6.1| ), the three- leg zig-zag ladder with PBC of type B can be considered as the 



natural counterpart of the N = 2 zig-zag ladder where the Hamiltonian can be written in 



the form ( 4.1 



In contrast to the preceding discussion of the two-leg case we will from now on again 
denote the length of each chain by L, i.e. the total number of spins is then 3L. Accordingly, 
we introduce a momentum p by a one-site translation along one of the three chains. 

In the following we will study each of these four boundary conditions in turn. 

A. The saturation field 

First, we compute the field h uc at which the transition to saturation takes place. This is 
not only useful to check the numerical computation to be reported in the next subsection, 
but like in the case N = 2 also serves as a guide where groundstates with incommensurate 
momenta can participate in the magnetization process. 

I) PBC of type A: 

A simple computation yields the excitation energy of one flipped spin above a ferromag- 
netic background as 
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£i«,±(p) = -2AJ' - — (l + cos(jj) ± sin(p)V3) - A J + Jcos(p) . (6.2) 

The minimum of this dispersion is located at 

tan min .,±) = ± j,_ 2J , (6-3) 

which is generically incommensurate. This then leads to the transition field 

J' 



h uc = -£is,±Gw,±) = A (2 J' + J) + j + V J' 2 ~ J' J + J 2 ■ (6.4) 

II) Open boundary conditions: 

The three-leg zig-zag ladder with open boundary conditions is shown in Fig. |l|. Here we 
consider only the case without dimerization (5 = 0). 

The dispersion of the gap corresponding to a single flipped spin above the ferromagnetic 
background is given by 
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f ls (p) = -AJ + Jcos(p)-J / (^A + ijA2 + 8cos(0 2 ) . (6.5) 



For AJ < J' < JVA 2 + 8, this dispersion has an incommensurate minimum at 



Pmin \ \ ^2.V 2 -2^.P 
cos — - — = . (6.6) 



2/4 J 

This then leads to an upper critical field 

/A \ 2 J' 2 3 J' A 

III) PBC of type B and C: 

Using the representation ( |6.1| ) one easily finds the dispersion of a single flipped spin 
above the ferromagnetic background 

£ u (p) = -A(2J' + J) + Jcos(p) + J' jcos (jt\ + cos (£j\ . (6.8) 

Since the PBC of type B (r = 2) are simpler to discuss analytically, we will now concentrate 
on this case. The dispersion (|6.8|) then always has an incommensurate minimum at 



/PminA ^JJ ' 2 " 3 JJ' + 9 J 2 - J' 

cos ( -^-j = . (6.9) 



From this we find an upper critical field for PBC of type B 



J' [9 J J' - 2 J' 2 + 27 J 2 ) ( J' 2 - 3 J J' + 9J 2 ^ 



K c = -S ls ( Pmhl .) = A (2J' + J) + — ^ — '- + ^ — '— . (6.10) 



B. Magnetization curves obtained by exact diagonalization 

Now we present one example of a magnetization curve for each of the four boundary 
conditions introduced above. In all cases we consider the SU(2) symmetric situation A = 1 
and set J' /J = 3/2. 

I) PBC of type A: 

An example of a magnetization curve for N = 3 with PBC of type A is shown in 
Fig. H For J = J', this geometry is one of the possible approximations to a two-dimensional 



triangular lattice antiferromagnet and a study of precisely this geometry |60| , |6ll exhibits 
a clear plateau with (M) = 1/3. For J = 0, one recovers an 'ordinary' rectangular spin 
ladder with periodic boundary conditions which then has equal coupling constants. The 
latter has already been investigated in some detail |IO| , p!3| , p!7| and evidence was found for a 
small plateau with (M) = 1/3 (there is further the possibility of a tiny spin-gap at h = 
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In Fig. |5| we chose J' = 3J/2 in order to study a new example. One can see from 
(3.3) that this gives rise to an incommensurability for (M) — > 1 (also at other values of the 



magnetization there is evidence that this choice of parameters lies in an incommensurate 
phase). Both at J = and J' = 0, one should choose even L to avoid frustration by the 
periodic boundary conditions along the chains, while for J' = J frustration is avoided if L 
is chosen to be a multiple of 3. To avoid both effects, L should be chosen as a multiple of 
6. However, this would considerably restrict the system sizes accessible to us. We therefore 
require only divisibility by three (which seems to be the more important one if the coupling 
constants are of a similar magnitude). 

The upper critical field for the choice of parameters in Fig. ^ is evaluated from ( |6.4j ) as 
h nc /J = 19/4 + y/1/2 « 6.0729. 

Like for N = 2 in Fig. [|, we applied the extrapolation procedure of |55| to the largest 
available system sizes in order to obtain the estimate for the magnetization curve in the 
thermodynamic limit which is shown by the bold line in Fig. [5]. 

In this extrapolated magnetization curve we have only drawn an (M) = 1/3 plateau but 
no further ones. This is motivated by the criterion ( |4.2| ): If we compare E(6) ~ 0.793 J with 
E(4) « 1.274 J (the latter is not shown in Fig. [5]) we see that the decrease is faster than 
( |4.2| ). These are only two data points which probably do not lie in the asymptotic regime (in 
particular non-monotonic finite-size corrections may still be important). It should also be 
noted that the weak-coupling analysis of Section II. A predicts a spin-gap. However, Fig. || 
corresponds to a rather large value J'. In any case, according to our numerical data, a 
spin-gap seems to be at least very small if present at all. We have therefore not drawn an 
(M) = plateau (corresponding to a spin-gap) in Fig. ||, though the speculative nature of 
the extrapolation in the region of small fields should be kept in mind. 

For non-zero magnetizations, we use the following generalization of ( |4.2j ) as the criterion 
for a vanishing plateau width (compare also ||63|| ) : 



h C2 (L)-h cl (L)~j-. (6.11) 

At (M) = 2/3 we find h C2 - h Cl « 0.993J, 0.584J, 0.529J and 0.331J for L = 4, 6, 8 and 12, 
respectively. Although we cannot entirely rule out a tiny (M) = 2/3 plateau on the basis 
of this data, it appears to be quite unlikely. We have therefore neither drawn a plateau at 
(M) = 2/3 in Fig. |5|. In any case, our main result for PBC of type A is that one can observe 
a clear plateau with (M) = 1/3. 

II) Open boundary conditions: 

The limit J — > of an iV-leg zig-zag ladder with open boundary conditions gives rise to 



the recently introduced 'diagonal ladders' [p4[ . In particular, for N = 3 with open boundary 



conditions one recovers the necklace ladder at J = 0. This necklace ladder is very similar 
to the S = (1, \) ferrimagnetic chain which is known to exhibit a plateau with (M) = 1/3 

US 



This strong-coupling limit J — > suggests to choose L divisible by two, as does the 
weak-coupling limit J' — > 0. On the other hand, at J' = J the N = 3 zig-zag ladder with 
open boundary conditions can be again considered as a strip of a triangular lattice which 
would suggest that L should be chosen a multiple of three. For the choice of parameters 
J' J J = 3/2 shown in Fig. |6] divisibility by three turns out to be not important and we 
impose only divisibility by two. 
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FIG. 5. Magnetization curves of an N = 3 zig-zag ladder with PBC of type A at J' = 3J/2. 
The lines are for L = 12 (full), L = 9 (dashed). L = 6 (dotted) and L = 3 (dashed-dotted). The 
bold line is an extrapolation to the thermodynamic limit. 



The choice of parameters corresponding to Fig. (A = 1, J' — 3 J/2) appears to lie 
entirely in an incommensurate phase. At least at the transition (M) — > 1, the momentum 
given by ( |6.6| ) is clearly incommensurate. The upper critical field is then found from fl6.7| ) 
as h uc = %J = 5.0625J. 

The finite-size magnetization curves in Fig. § exhibit the expected plateau at (M) = 1/3. 
In addition, there is clear evidence for a further plateau at (M) = 2/3. The boundaries of 
these two plateaux were extrapolated applying either a Shanks transform (which is the a = 
special case of the vanden Broeck-Schwartz algorithm - see e.g. p7| ) to their values h Ci (L) 
at the available system sizes L, or (in the case of the upper boundary of the (M) = 1/3 
plateau) by fitting to 



h Ci (L) = /i Cl (oo) 



a 

L 



(6.12) 



This heuristic formula is motivated by ( |6.11| ), since eq. ( |6.12| ) can be expected to yield 
coincident plateau boundaries if there is actually no plateau. In general, (|6.12|) will work well 
if it is applied to system sizes substantially below the correlation length and underestimate 
the width of a plateau otherwise. Here, it gives a small correction to the value for h c . obtained 
for the largest L, as does the Shanks extrapolation at the other plateau boundaries. 
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FIG. 6. Magnetization curves of an iV = 3 zig-zag ladder with open boundary conditions 
along the rungs and J' = 3 J/2. The lines are for L = 12 (full), L = 10 (dashed). L = 8 (dotted), 
L = 6 (long dashed-dotted) and L = 4 (dashed-dotted). The bold line is an extrapolation to the 
thermodynamic limit. 



The finite-size spin-gap for L = 4, 6 and 8 is nicely fitted by ( |4.2|) . This indicates that 
there is no spin-gap at h — and we have therefore not drawn an (M) = plateau in the 
extrapolated magnetization curve of Fig. |[ Recently it was shown that a related three-leg 
ladder at zero field is also massless and gives rise to a c = 2 (su(2) 2 x Ising) conformal 



field theory |68 |. However there the chirally asymmetric perturbation was eliminated by a 
particular choice of coupling in order to permit an analytic treatment. In the present case 
the central charge might therefore be smaller than two. 

The bold line in Fig. |6| is an extrapolation taking into account the foregoing discussion 
of plateau boundaries. Between the plateaux it has been obtained by the same procedure 
as used in Figs. § and |]. 

Ill) PBC of type B and C: 

Inspection of the r = 2 version of ( |6.1|) shows that at J = 0, the iV = 3- leg zig-zag ladder 
with PBC of type B is equivalent to the N = 2-leg ladder with equal coupling constants. So, 
the results for the two-leg case discussed earlier can be carried over to this three-leg ladder 
at strong coupling. For example, this iV = 3-leg ladder should exhibit incommensurate 
groundstate momenta at J' ^> J. More important for our purposes is that a small gap, but 
no other non-trivial plateaux are expected at strong coupling. To see to which extent this 
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is generic, we show in Fig. [j] numerical results obtained for the smaller value of J 1 = 3J/2. 
Fig. H shows the analogous result for PBC of type C. 

Since at J = J', neither PBC of type B nor PBC of type C respect the sublattice 
structure of the triangular lattice, there is no reason to expect those L which are a multiple 
of three to play any particular role. We therefore simply consider even L. 
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FIG. 7. Same as Fig. |, but for PBC of type B. 



If one ignores the bold extrapolated curves in Figs. [7| and |8|, the finite-size magnetization 
curves look quite similar at least if compared to the two boundary conditions discussed 
before. Such a similarity of PBC of type B and type C is also suggested by the fact that 
both of them can be mapped to a single chain QS.1| ). However, there are definitely at least 
quantitative differences, as already the value of the upper critical field shows: At A = 1, 
J' J J = 3/2, one finds from (|6~Tol) that h uc /J = 5 + §v% « 5.6495 for PBC of type B, while 
minimization of ( |6.8|) with r = 4 yields h uc m 6.1 193 J for PBC of type C. 

The finite-size data for L = 4, 6 and 8 indicates a small spin-gap at zero field. For 
PBC of type B, application of a 1/L-extrapolation in the spirit of ( 6.12J) leads to a gap 
E « 0.24J while the Shanks transform yields E « 0.64J. Though both methods agree 
on the presence of a gap, we are not able to determine it accurately - by comparison we 
conclude E/J = 0.44 ± 0.20 for PBC of type B. Similarly, for PBC of type C, the Shanks 
transform yields E « 0.27J while from ( |6.12| ) one finds E rj 0.48J. Here, we use the latter 
value though a large uncertainty should be kept in mind. 

While one observes clear plateaux with (M) ^ in Figs. ||] and [| it is not immediately 
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FIG. 8. Same as Fig. §, but for PBC of type C. 



clear from the finite-size data of Figs. [?| and |8] if such non-trivial plateaux survive the 
thermodynamic limit. This issue therefore requires a more detailed discussion. First we 
take a closer look at (M) = 1/3. For PBC of type B, the finite-size plateau with L = 4, 6 
and 8 is roughly (but not very well) fitted by ( |6.11|) . More strongly, one finds a negative 
plateau width if one tries to apply the extrapolation formula ( |6.12| ) to its boundaries. This 
indicates that there is not even a tiny plateau with (M) = 1/3 in the thermodynamic limit 
L — > oo. For PBC of type C in contrast, the plateau width at L = 8 is actually larger than 
that at L — 6. We have therefore drawn an (M) = 1/3 plateau in Fig. |8] whose boundaries 
have been determined using (|6.12|) . 

At (M) = 2/3 we have data for a further system size [L = 10). If we discard the L = 4 
data and ignore some apparent non-monotonic finite-size effects we actually obtain fairly 
good agreement with ( |6 . 1 1| ) in both cases. So, we do not find a plateau at (M) = 2/3 in 
either of the two cases. 

Extrapolated magnetization curves are shown by the bold lines in Figs. and §. They 
have some remaining wiggly features which are related to the fact that non-monotonic finite- 
size effects may still be important at the system sizes used for the extrapolation. 

The absence of a plateau e.g. at (M) = 1/3 can be attributed to enhanced translational 
symmetry: The N = 3 zig-zag ladder with PBC of type B can be mapped to a single 
Heisenberg chain (|6.1|) . A plateau at (M) = 1/3 would then require breaking of the enhanced 
translational symmetry by / = 3 sites and a plateau at (M) = 2/3 would require the even 
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larger period of / = 6 sites (compare ( |1 . 1|) where V = I due to the enhanced symmetry). 
However, to the best of our knowledge spontaneous symmetry breaking with periods of more 
than two sites has not been observed so far (at least for ordinary ladders it is unlikely to 
occur fll3| ). 

It should be mentioned though that the same argument can be applied to PBC of type 
C, where a plateau with (M ) = 1/3 seems likely Furthermore, all versions of PBC are very 
similar from the weak-coupling point of view (compare Section III). In particular, one would 
expect a plateau with (M) = 1/3 for all versions of PBC if it appears for one of them in the 
weak-coupling region. Which plateaux are actually absent and which ones present therefore 
needs further investigation. 

VII. DISCUSSION AND CONCLUSION 

In this paper we have analyzed zig-zag coupled chains using a range of field-theoretical 
and numerical techniques. First, we have discussed the case of zero field and weak inter- 
chain coupling where in a field-theoretical formulation the most relevant interaction for the 



ordinary ladders is replaced with a chirally asymmetric one [24 1. In the SU(2) symmetric 
case this chirally asymmetric term is marginal and then current- current interaction terms 
(usually not considered in the presence of the relevant terms that arise in the normal ladders) 
have to be included in the RG analysis, leading us to our one-loop RG equations fl2.10|) . In 
this approximation, one sees that the effect of the chirally asymmetric term is to push the 
system further into the massive phase for the AF case, and to open a gap for small J' < 0. 
The same is true for N > 3, when all the chains are (weakly) antiferromagnetically coupled 
in an equal manner. For sufficiently small A < 1, the zero magnetization groundstate is a 
massless c = 1 theory, in contrast to the N = 2 ladder with normal couplings. 

We have then shown that this more relevant interaction is restored by a magnetic field, 
dimerization or doping with charge carriers. The latter variants of zig-zag ladders should 
therefore be similar to the usual spin ladders. We have also analyzed the appearance of 
plateaux in the magnetization curves of N > 3 zig-zag ladders. We have found that the 
situation is similar to that encountered in the normal ladders (in the weak coupling limit) 
|T3|j , except for a trivial rescaling of the couplings which lead to slightly different quantitative 



predictions for the opening points. 

These considerations may also be relevant to the description of experiments. Firstly, we 
have seen that under very general conditions the normal perpendicular couplings are those 
that really matter. Secondly, we have argued that dimerization along the legs (see Fig. 
0) might be an important feature in understanding the magnetization experiments JZS 



on 



NH4CUCI3 (see also p(J for a related observation in the context of the usual spin ladders 
with dimerization along the chains). In the present context, this is suggested by the limit 
of strong dimerization <5 — >- 1 which we have discussed in the appendix. Dimerization also 
seems to be a small enough modification to be plausible, but this will have to be confirmed 
by a determination of the structure of NH4CUCI3 at low temperatures. 

The field-theoretical approach is complemented in a second part by a numerical analysis. 
First, we have analyzed a two- leg zig-zag ladder without dimerization. For positive J' 
we have confirmed known facts such as the appearance of a spin-gap but no non-trivial 
plateaux in the magnetization curve. The magnetized groundstates exhibit two different 
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types of behaviour: At small J' incommensurate groundstate momenta participate in the 
magnetization process, while for small J all momenta are commensurate. For the case of 
ferromagnetic coupling J' < we have found interesting new behaviour: At small \J'\ spins 
flip in pairs and for this reason all magnetized groundstates are commensurate. At larger 
| J' | more than two spins can flip simultaneously and at the same time also incommensurate 
groundstate momenta become relevant to the magnetization process. Finally, for J 7 < — 4 J, 
the behaviour is completely ferromagnetic. The 'phases' with incommensurate groundstate 
momenta are reminiscent of similar phenomena observed in a spin-one chain |50| , |5T[| . It has 



been conjectured [ ]5U| , |o"3|] that also the two-leg zig-zag ladder gives rise to a two-component 
Luttinger liquid, but at least in the weak-coupling region we did not find evidence for more 
than one massless degree of freedom. 

We have then investigated the transition to saturation in more detail and confirmed the 



univeral DN-PT behaviour (|5.5|) . The only exception is J' = 4 J where a quartic behaviour 
( p.6| ) was already observed in [f4Tj| . This modified exponent has a simple explanation in 
terms of an exceptional behaviour of the single-spinwave dispersion which appears not to 
have been noted before. In the region of not too strong ferromagnetic J' < spins flip in 
pairs and this relation to the single-spinwave dispersion is lost. It is therefore non-trivial 
that our numerical data for the asymptotics of the magnetization curve is still consistent 
with the DN-PT universal behaviour. 

Finally, we have computed magnetization curves for four different variants of a three-leg 
ladder. We observed a remarkable dependence on the geometry of the interchain coupling 
(see Figs. §-§). Clear plateaux can be observed in several cases while indications against 
the presence of such plateaux were obtained in other situations []. We expect that the same 
types of commensurate and incommensurate phases as observed for the two-leg ladder are 
also present in zig-zag ladders with three or more legs. This is indicated e.g. by the analysis 
of single-spinwave excitations for three-leg ladders, but would need verification. Further new 
phenomena might appear for more than two coupled chains. 

All the observed plateaux can be interpreted in terms of the quantization condition ( |1 . 1| ) . 
The essence is enhanced translational symmetry in certain cases which induced by frustration 
is then spontaneously broken (in most cases to periods I < 2). Specifically, for both the 
ordinary and the zig-zag two-leg ladder at 5 = one then has a translationally invariant unit 
cell containing V = 2 spins. However, the interpretation is different: For zig-zag coupling, 
translational symmetry is first enhanced by a factor of two and then spontaneosly broken by 
a period I = 2. This also explains why one needs a dimerized interchain coupling (0 < 5 < 1) 
to have a plateau with (M) = 1/2 in the two-leg zig-zag ladder [|l|,[^,f|(| which is permitted 
by ( |1 . 1[ ) with N = 2, I = 2 only in the presence of dimerization. Dimerization along the 
legs breaks translational symmetry further to V = 41. Frustration- induced spontaneous 
symmetry breaking to a period I = 2 then permits the aforementioned appearance of further 
plateaux with (M) = 1/4 and 3/4. 

To interpret our results for the three- leg zig-zag ladders in terms of ( |1 . 1|) , one should 



3 It should be noted that due to non-monotonic finite-size effects it is difficult to reliably exclude 
plateaux not only in the cases discussed in Section VI but also otherwise. Nevertheless, the evidence 
for or against a plateau e.g. at (M) = 1/3 for N = 3 depends strongly on the boundary conditions. 
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substitute V = 31 for OBC and PBC of type A and V = I for PBC of type B and C. With 
one exception, all the plateaux sketched in Figs. |5|-|8| can then be naturally interpreted with 
periods I < 2. Only the plateau with (M) = 1/3 for PBC of type C in Fig, |8| requires 
a period I = 3 which (if the presence of this plateau is confirmed) would be the highest 
observed period which we are aware of. 

In fact, the numerical support for a plateau with (M) = 1/3 for N = 3 and PBC of type 
C is rather weak and would deserve further attention. Conversely, the absence of a spin- 
gap in Fig. |^ is not on safe grounds, and the Abelian bosonization analysis would actually 
predict at least a small spin-gap in the weak-coupling region. Another issue for N = 2 
which requires further attention is that there is no numerical evidence yet for the spin-gap 
predicted by field theory for weak ferromagnetic J' < 0. 

So far, we have just observed numerically for N = 3 that modifications in the boundary 
conditions have a drastic effect on the magnetization process. This is in contrast to ordinary 
spin ladders [|H| and not apparent in the field theoretical treatment. It may be necessary to 
include higher loop corrections or to perform a non-perturbative analysis of the RG equations 
in order to understand how the novel interaction term arising from the zig-zag coupling gives 
for example rise to the differences in the different versions of PBC. 

In summary, on the one hand we believe that we have exhibited interesting properties of 
zig-zag spin ladders. On the other hand, there is a number of points which deserve further 
attention. We hope that with this combination, the present paper will stimulate further 
research on zig-zag spin ladders. 
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APPENDIX A: STRONG-COUPLING EFFECTIVE HAMILTONIAN FOR A 

DIMERIZED TWO-LEG LADDER 

In this appendix we consider the Hamiltonian for a two-leg zig-zag ladder with dimerized 
chains and coupling between the chains: 

i=l x=l 

+f E {(! + <*')<V ■ 4* + (i - <*')&,* • &,*+i} ( A1 ) 

x=l 
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If we assume that J'(l ± 5'),J(1 — 5) <C J(l + 5), we can describe the transition from 
(M) = to (M) = 1 by an effective Hamiltonian following 0,I2,|TB-|T^|. On each 'rung' (a 



bond coupled with coefficient J(l + 6)) we retain only two states: The singlet and the fully 
polarized state. With an appropriate choice of basis, the first-order effective Hamiltonian 
can then be written as (note that we include the external magnetic field h in zeroth order): 



( 1 1 

-TCeff. = ^(1 — ") Z^ j T \px ^x+2 + b x O x+2 ) + ~7b x b x -\ 

+ f (1 " f) £ {; (S}A +1 + ^S+ +1 ) + j^At,} 



+ f (1 + 3*') £ i (5+ +1 5 2 " +2 + S 2 " +1 S+ +2 ) + J' (3 + <?) £ is 2r+1 S 



+(^-*>4)(i+|:4 



(A2) 



The first line comes from one weak coupling J (1 — 5) along the two original chains, the 
second one from one coupling J' (1 — 5') between the two chains and the third line arises 
from one coupling J' (1 — 5') plus two couplings J' (1 + 5') between the chains. The fourth 
line is just a first-order correction to the external magnetic field (apart from a trivial additive 
constant). 

In eq. (|A2|) one recognizes again a two-leg zig-zag ladder with dimerized coupling between 
the chains. The coupling along the two chains is J(l — 5) with an effective XX Z anisotropy 
A e ff. = 1/2. The coupling between the two chains is not only dimerized but also has an 
alternating XXZ anisotropy: J' cS = J'{\ — 5'), A' off = 1/2 on even sites, J' eS = J'{\ + 35'), 
^off. = (3 + 5')/ (2(1 + 35')) on the odd sites. For these values of parameters, the two- 
leg zig-zag ladder in a magnetic field has unfortunately not yet been studied in detail. 
However, this mapping is still suggestive since we know that a two-leg zig-zag ladder with 
dimerized coupling between the chains can exhibit plateaux at (M e ff.) = and (M e g.) = ±1/2 
r2L|3^flIJ . This suggests the possibility of the two-leg ladder with dimerized coupling along 



the chains (|A1|) having not only plateaux at (M) = and (M) = 1/2 (corresponding to 
(M eff .) = 0) but also at (M) = 1/4 ((M eff .) = -1/2) and (M) = 3/4 ((M eff .) = 1/2). This 
observation may be relevant to the magnetization experiments on NH4CUCI3 f28f , since these 
are precisely the observed magnetization plateaux. However, in the strong dimerization limit 
considered in this appendix, one would certainly have a pronounced (M) = plateau which 
is not observed in NH4CUCI3. So, even if dimerization along the chains should be important 
in this compound, it cannot really lie in the region where the above mapping is applicable. 
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